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ABSTRACT 

Several interesting correlations among Gamma Ray Burst (GRB) observables 
with available redshifts have been recently identified. Proper evaluation and 
calibration of these correlations may facilitate the use of GRBs as standard 
candles constraining the expansion history of the universe up to redshifts of 
z > 6. Here we use the 69 GRB dataset recently compiled by Schaefer (2007) 
and we test the calibration of five of the above correlations (1 : Ep^ak — Ej, 
2 : Epeak — L, 3 : Tiag — L, A : V — L, 5 : trt — L) with respect to two potential 
sources of systematics: Evolution with redshift and cosmological model used in 
the calibration. In examining the model dependence we assume flat ACDM and 
vary Q.^. Our approach avoids the circularity problem of previous studies since 
we do not fix to find the correlation parameters. Instead we simultaneously 
minimize with respect to both the log-linear correlation parameters a, b and 
the cosmological parameter fim. We find no statistically significant evidence 
for redshift dependence of a and b in any of the correlation relations tested [the 
slopes of a(z) and b{z) are consistent with at the 2(t level]. We also find that 
one of the five correlation relations tested {Epeak — E^) has a significantly lower 
intrinsic dispersion compared to the other correlations. For this correlation 
relation, the maximum likelihood method (x^ minimization) leads to 6i = 
50.58 ± 0.04, ai = 1.56 ± 0.11, = O.SOlJ;^^ respectively. The other four 
correlation relations minimize for a flat matter dominated universe flm — 
1. Finally, a cross-correlation analysis between the GRBs and Snia data for 
various values of fim has shown that the Ep^ak — E-y relation traces well the 
SnIa regime (within 0.17 < z < 1.755). In particular, for flj^ — 0.15 and 
rim — 0.30 we get the highest correlation signal between the two populations. 
However, due to the large error-bars in the cross-correlation analysis (small 
number statistics) even the tightest correlation relation {Ep^ak — E-y) provides 
much weaker constraints on than current SnIa data. 
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1 INTRODUCTION 

Several cosmological observations (Riess et al. 1998; 
Perlmutter et al. 1999; Tegmark et al. 2006; Davis et 
al. 2007, Komatsu et al. 2008) have converged dur- 
ing the last decade towards a cosmic expansion his- 
tory that involves a recent accelerating expansion of 
the universe. This expansion has been attributed to an 
energy component (dark energy) with negative pres- 
sure which dominates the universe at late times and 
causes the observed accelerating expansion. The sim- 
plest type of dark energy corresponds to the cosmologi- 
cal constant (Padmanabhan 2003). Alternatively, mod- 
ifications of general relativity have also been utilized 
to explain the observed cosmic acceleration (Boisseau 



et al. 2000; Perivolaropoulos 2005; Caldwell, Cooray & 
Melchiorri 2007; Nesseris & Perivolaropoulos 2006; Jain 
& Zhang 2007; Wang et al. 2007; Heavens, Kitching & 
Verde 2007; Nesseris & Perivolaropoulos 2007). 

The geometrical probes (Lazkoz, Nesseris & 
Perivolaropoulos 2008) used to map the cosmic expan- 
sion history involve a combination of standard can- 
dles [Type la supernovae (SnIa) Davis et al. 2007] and 
standard rulers [clusters, CMB sound horizon detected 
through Baryon Acoustic Oscillations (BAO; Percival 
et al. 2007 ) and through the CMB perturbations an- 
gular power spectrum Komatsu et al. 2008]. These ob- 
servations probe the integral of the Hubble expansion 
rate H{z) either up to redshifts of order 2; ~ 1 — 2 
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(Snia, BAO, clusters) or up to the redshift of recom- 
bination {z ~ 1089). Alternatively, dynamical probes 
(Bertschinger 2006; Ncssoris & Perivolaropoulos 2008) 
of the expansion history based on measures of the 
growth rate of cosmological perturbations arc also con- 
fined to relatively low redshifts up to z ~ 1. It is there- 
fore clear that the redshift range ~ 2 — 1000 is not di- 
rectly probed by amy of the above observations. Even 
though many models of dark energy predict a decel- 
erating expansion in that redshift range due to mat- 
ter domination, the possibility of non-trivial expansion 
properties at higher redshifts can not be excluded. In 
order to investigate this possibility we need a visible 
distance indicator at redshifts z > 2. 

Gamma Ray Bursts (GRBs) are the most violent 
and bright explosions in the universe. They are pro- 
duced by a highly relativistic bipolar jet outflow from a 
compact source (Rhoads 1999; Piran 2004). At present 
the most distant GRB (GRB 050904) is at a redshift 
z = 6.3 (Kawai et al. 2005). The fact that GRBs are 
detected up to very high redshifts makes it tempting 
to try and use them as standard candles that could be 
used to constrain the cosmological expansion history in 
a similar way as SnIa. The problem is that GRBs ap- 
pear to be anything but standard candles: they have 
a very wide range of isotropic equivalent luminosities 
and energy outputs. Several suggestions have been made 
to calibrate them as better standard candles by using 
correlations between various properties of the prompt 
emission, and in some cases also the afterglow emission. 
While there is good motivation for such cosmological ap- 
plications of GRBs, there are many practical difficulties. 
Indeed, a serious problem that hampers such a straight 
forward approach is the intrinsic faintness of the nearby 
events which introduces a bias towards low (or high) 
values of GRB observables and as a result of this, the 
extrapolation of each of the correlations to low-z events 
is faced with serious problems. An additional problem 
is that the GRB surveys suffer, due to the unknown flux 
limit, from the well known degradation of sampling as a 
function of redshift (Lloyd, Petrosian, & Mallozzi 2000). 
One might also expect a significant evolution of the in- 
trinsic properties of GRBs with redshift (also between 
intermediate and high redshifts) which can be hard to 
disentangle from cosmological efi^ects. In addition, even 
after properly accounting for the observed correlations, 
the scatter in the luminosity of the standardized can- 
dles is still fairly large. Finally, the calibration of the 
observed correlations require the assumption of a cos- 
mological model (luminosity distance vs redshift) in or- 
der to convert the observed bolo-mctric peak flux or 
bolo-metric fluence to isotropic absolute luminosity L 
or to a total coUimation corrected energy Ej. The use 
of a cosmological model to perform the calibration cre- 
ates a circularity problem and a model dependence of 
the obtained calibration. 

Despite of the above difficulties, the potential ben- 
efits of obtaining even approximate standard candles at 
redshifts as high as 2 > 6 has prompted a siginficant 
activity towards both testing the usefulness of GRBs as 
standard candles (Amati et al. 2002; Ghirlanda, Ghis- 
ellini & Firmani 2006; Firmani et al. 2006; Li 2007a; Li 



2007b; Butler et al. 2007; Zhang & Xie 2007; Hooper & 
Dodelson 2007) and eagerly utilizing them to constrain 
cosmological parameters (Schaefer et al. 2003; Zhang & 
Meszaros 2004; Dai, Liang & Xu 2004; Di Girolamo, et 
al. 2005; Schafer 2007; Bertolami & Tavares Silva 2006; 
Wang & Dai 2006; Demianski, et al. (2006); Li et al. 
2008a; Amati et al. 2008) . This activity has lead to a de- 
bate about the usefulness of GRBs as standard candles 
with both discouraging (Li 2007a; Li 2007b; Butler et 
al. 2007) and encouraging (Amati et al. 2002; Ghirlanda 
et al. 2006; Firmani et al. 2006; Schaefer 2007; Zhang & 
Xie 2007; Hooper & Dodelson 2007) results. 

The goal of the present study is to evaluate the cur- 
rent utility of GRBs as standard candles and identify the 
correlation relations that are more promising in the de- 
termination of the isotropic absolute luminosity and to- 
tal coUimation corrected energy of the GRBs. The eval- 
uated correlation relations involve the relationship be- 
tween a measurable observable of the light curve (lumi- 
nosity indicator) with the GRB luminosity, given by the 
correlation relation in the form of a power-law (Schaefer 
2007), ic., E., = B^, peak Epeak"-' (Ghirlanda, Ghis- 
elhni & Lazzati 2004) L — BpeakEpeak"'"'^'' , (Schae- 
fer et al. 2003) L = BiagTiag"'"'' , (Norris, Marani & 
Bonnell 2000) L = B^V^ , (Fenimore & Ramirez-Ruiz 
2000) and L = Brttrt"^'^ (Schaefer 2007) [see section 
2 for the definitions of the above observable luminosity 
indicators]. These five correlation relations can be de- 
noted compactly a.s R = Bj Q'^i where R is proportional 
to the GRB absolute luminosity L, Q is a GRB observ- 
able, j = 1, .., 5 counts correlation relations while Bj, aj 
are parameters to be calibrated using a minimization 
with the GRB data. 

We use the 69 GRB dataset compiled by Schae- 
fer (2007) and fit the logarithmic linear form of the 
correlation relations ie., log_R = bj + aj log Q (where 
bj = log Bj ) using a maximum likelihood method of 
minimization (symmetric in R and Q). We implement 
the following tests: 

• Wc compare the quality of fit obtained with each 
one of the above five correlation relations by evaluating 
Xmin per degree of freedom in each case. 

• We split the dataset into four consecutive redshift 
bins and evaluate the best fit parameters aj , bj in each 
bin to test for evolution effects of each correlation rela- 
tion. 

• The value of R (e.g. isotropic absolute luminosity) 

depends on both the observed GRB bolo-metric peak 
flux and on the cosmological model used to evaluate 
the luminosity distance (see eq. 1). We thus evaluate 
the dependence of the quality of flt {Xmin) on the as- 
sumed cosmic expansion history H{z). For concreteness 
we assume a flat ACDM model and fit aj, bj for various 
Jim. This test is equivalent to simultaneously minimiz- 
ing with respect to both the calibration parameters 
aj and bj and the cosmological parameter Om. Such an 
approach avoids the circularity problem discussed above 
and leads to an unambiguous determination of the cos- 
mological parameter Sim. We would like to stress here 
that Schaefer (2007) claims that during the fim fitting 
a marginalization over the parameters aj, bj could po- 
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tentially solve the circularity problem. However, if the 
marginalization is done over the la range of the best fit 
parameter values then the information for the imposed 
value of fim is still carried over leading to circularity. If 
it is done on a wider range then the errors on fJm would 
be even larger than the ones obtained here. 

• Finally, we perform a direct cross-correlation com- 
parison of the GRB data Schaefer (2007) with a recent 
Snia dataset Davis et al. (2007) in the redshift range of 
overlap. 

The structure of this paper is the following: In the 
next section we briefly review the five correlation rela- 
tions studied and describe our numerical fitting meth- 
ods. In section 3 we present our numerical results with 
respect to both the redshift dependence of the correla- 
tion parameters and their dependence on the assumed 
cosmological model. In section 4 we compare the GRBs 
with the SnIa data in the redshift range of overlap. 
Finally, in section 5 we conclude and discuss future 
prospects of this work. 



2 CORRELATION RELATIONS - METHOD 

The correlation relations discussed in what follows con- 
nect GRB observables with the isotropic absolute lumi- 
nosity L or the coUimation corrected energy E~, of the 
GRB. Such observable properties of the GRBs include 
the peak energy, denoted by Epeak, which is the pho- 
ton energy at which the i' F,^ spectrum is brightest; the 
jet opening angle, denoted by d-^c^t, which is the rest- 
frame time of the achromatic break in the light curve 
of an afterglow; the time lag, denoted by Tiag, which 
measures the time offset between high and low energy 
GRB photons arriving on Earth and the variability, de- 
noted by V, which is the measurement of the "spikiness" 
or "smoothness" of the GRB light curve. In the litera- 
ture, there is a wide variety of choices for the definition 
of V (Fenimore & Ramirez-Ruiz 2000; Reichart et al. 
2001). In this work, wo follow the notations of Schae- 
fer (2007) in which the observed V value varies as the 
inverse of the time stretching, so the corresponding mea- 
sured value must be multiplied by 1 + z to correct to the 
GRB rest frame. An additional luminosity indicator is 
the minimum rise time Schaefer (2007) denoted by trt, 
and taken to be the shortest time over which the light 
curve rises by half the peak flux of the pulse. 

These quantities appear to correlate with the GRB 
isotropic luminosity or its total coUimation-corrected 
energy. This property can not be measured directly but 
rather it can be obtained through the knowledge of ei- 
ther the bolo-metric peak flux, denoted by Pboio; or 
bolo- metric fluence; denoted by Sboio, Schaefer (2007). 
Therefore, the isotropic luminosity is given by: 

L = 4ndl{flm,z)Pi,oio (1) 
and the total coUimation-corrected energy reads: 

= A-Kdliflm, «)S'boloFbeam(l + «)"^ (2) 



where Fbeam is the beaming factor* (1 — cosSjet). 
The luminosity correlation relations arc power-law rela- 
tions of either L or Ej as a function of nag, V, Epeak, 
Trt- Both L and E^y depend not only on the GRB ob- 
servables Pboio or Sboio, but also on the cosmological pa- 
rameters through the luminosity distance dL{z) which 
in a flat universe is expressed in terms of the Hubble 
expansion rate H{z) = HoE{z) as 



dL{Clm,z) = {1 + z) — 

no 



dz' 



where E'^(z) = f2m(l -|- z)^ + iixfx{z) and the dimen- 
sionless dark energy density /x(a) is given by 



Px(z) 

Px(0) 



(3) 



where the last equality is valid in the case of ACDM 
(fix = Oa) which is assumed in what follows. 

The relationship between a measurable observable 
of the light curve (luminosity indicator) with the GRB 
luminosity is given by the luminosity relation in the 
form of a power-law, i. e., E^ = Bj^peakEpeak""''''"'''' 
(Ghirlanda et al. 2004) L = BpeakEpeak'^"'"''' , (Schae- 
fer et al. 2003) L = Biagr/ag"'"' , (Norris et al. 2000) 
L = (Fenimore & Ramirez-Ruiz 2000) and 

L = Brttrt""'^ (Schaefer 2007). The observed (on 
Earth) luminosity indicators will have different values 
from those that would be estimated in the rest frame 
of the GRB event. That is, the light curves and spectra 
seen by the Earth-orbiting satellites suffer time dilation 
and red-shifting. Therefore, the physical connection be- 
tween the indicators and the luminosity in the GRB rest 
frame must take into account the observed indicators 
and correct them to the rest frame of the GRB. For 
the temporal indicators, the observed quantities must 
be divided by 1 -|- a to correct the time dilation. The 
observed V- value must be multiplied by 1 + 2: because it 
varies inversely with time, and the observed Epeak must 
be multiplied by 1 -I- z to correct the redshift dilation of 
the spectrum. We have also rescaled the luminosity in- 
dicators to dimensionless quantities by using the same 
values employed by Schaefer (2007) in order to mini- 
mize correlations between the normalization constant 
and the exponent during the fitting, i. e., for the tem- 
poral luminosity we use 0.1 second, for the variability 
0.02, and for the energy indicator SOOkeV. For example 
the effective GRB frame dimensionless value of Ep^ak 
used in our analysis is Q = ^^^Q^^^.y'^^ where Ep^ak is 
obtained from Table 1. The 69 GRB data of Schaefer 
(2007) used in our analysis arc shown in Table 1. 

To explain the calibration procedure in general, we 
denote the five luminosity relations hy Ri = B and 



* Note, that -Fbeam is calculated with the aid of di^ (z) (see 
Sari, Piran & Halpcrn 1999; Butler, Kocevski, & Bloom 
2008). For small values of 0jet, the dependence of Fi,eam 



.,-1/2 



Using the 



on dij{z) is found to be: Fbeam "-j^ 
data from Schaefer 2007 (Oin = 0.27) we have to multi- 
ply the corresponding Ff,eam value by the following factor: 

[di(0.27, z)/dL{nrn,z)]^/''. 
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Table 1. The dataset of GRB observables used in our analysis (from Schaefer 2007). 
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GRB 



^'bolo 

[erg/cm^s] 



[erg/cm^ 



'''lag 

[sec] 



V 



-^pcak 

[keV] 



TRT 
[sec] 



060115 3.53 1.3E-7 ± 1.6E-8 

060116 6.60 2.0E-7 ± l.lE-7 

060124 2.30 l.lE-6 ± 1.2E-7 

060206 4.05 4.4E-7 ± 1.9E-8 

060210 3.91 5.5E-7 ± 2.2E-8 

060223 4.41 2.1E-7 ± 3.7E-8 

060418 1.49 1.5E-6 ± 5.9E-8 

060502 1.51 3.7E-7 ± 1.6E-7 

060510 4.90 l.OE-7 ± 1.7E-8 

060526 3.21 2.4E-7 ± 3.3E-8 

060604 2.68 9.0E-8 ± 1.6E-8 

060605 3.80 1.2E-7 ± 5.5E-8 

060607 3.08 2.7E-7 ± 8.1E-8 



3.37E-5 
1.94E-5 



± 3.4E-6 
± 1.2E-6 0, 



1.17E-6 ± 1.7E-7 0, 
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we take their logarithms to express them as a hnear 
relation of the form 



log J?i = log B + o log Qi ^ j/i = 6 + o Xi 



(4) 



where yi = log Ri, Xi = log Qi and b = logB. Notice 
that Ri depends on the cosmological parameters [eg. 
nin) through the luminosity distance (see eqs (1) and 
(2)1. 

In order to find the best fit values of the parameters 

a and b wc use a symmetrized form of the maximum 
likelihood method and minimize {a,b,^lui) defined as 



X^(o,6, fin 



[2/i(f2m) - ytia,b)f 



2 

(J sys ' 



(5) 



Note that each correlation point is weighted by its 
error" which means that points with large errors have a 
negligible contribution to the x^i^-j ''i ^m) function. Fur- 
thermore, we normalize over the number of points^ 
and since the number of points is not too low (see last 
column of Table 2) we do not anticipate the number of 
points by itself to introduce a further significant bias. 
The quantities appearing in eq. (5) are defined as fol- 
lows: 

• 2/i(f2m) corresponds to either logL or logS-^ as ob- 
tained from Table 1 and eqs. (l)-(2) and depends on 
the cosmological parameter fim through the luminosity 
distance dL{^m, Zi). 

• yi{a,b) is the predicted value of j/i = logi?, on the 
basis of the linear logarithmic relation (4). 

• Using error propagation we find (7y^ from the errors 
of Table 1 as 



dy{R^) 
dR, 



In 10 Pboio 



(6) 



when y corresponds to log L and 



t Notice tiiat the corresponding errors assumed to follow a 
normal distribution, which is the usual requirement in order 
to have a distribution. 



lnl0\/( 



\ ^ beam / 



(7) 



when y corresponds to log Ej. In the case of asymmetric 
errors we utilize the simplest approach of symmetrizing 
the errors by taking their mean value on the two sides. 
Similarly, we find as 



dx{Qi) 



dQi 



In 10 Qi 



(8) 



scatter determined for each correlation relation. Note 
that we do not treat aays as a free parameter (we do 
not minimize with respect to it) but wc fix it only at 
the end of the minimization by demanding Nx^ /dof ~ 
= 0{1). Thus we treat asys not as a physics related 
parameter but as an unknown source of scatter which 
is required to make the quality of fit acceptable. A sim- 
ilar approach has been used by the SNLS collaboration 
(Astier et al. 2006) in cosmological fits of type la super- 
nova data. This approach however should be generally 
used with care (Vishwakarma 2007) since it may hide a 
low quality of distance indicators or it may introduce a 
good quality of fit for a model, by brute force. 

Equation (5) has the advantage of being symmetric with 
respect to errors in both the Xi and yt variables (see eg 
Amati et al. 2008 and references therein). It is interest- 
ing to mention here that the statistical results depend 
on where we put the asys parameter (either in x, y-axis 
or both). For example asys could have been included 
along with the x-axis error as 



X^(o,6,fii„) = - 
^ — ' a. 



M^m) - yi{a,b)f 
+ a sys + cL^{(^li + a'iys) 



(9) 



leading to full symmetry between x and y. We have ver- 
ified immcrically that the use of eq. (9) [proposed by N. 
Butler private communication] leads to values of best 
fit parameters a and 6 that depend very sensitively on 
the value of asys- We consider this to be an undesir- 
able feature and thus we have chosen to use the more 
robust form of (5). Schaefer (2007) on the other hand, 
uses a linear regression procedure by putting asys in 




both X and y. The novelty of our statistical analysis is 
that we treat the problem with a maximum likelihood 
which includes the data errors and thus the correspond- 
ing results are less dependent on the assumption related 
with the value of the ffsys parameter. A symmetrized 
method may be warranted because the two variables Ri, 
Qi are not directly causative [for example, it is possible 
that the scatter may be dominated by hidden variables, 
or variables not directly measured or treated, such as 
the bulk Lorentz factor Fjet (Schaefer et al. 2003, Am- 



ati et al. 2002)]. The bisector method of two ordinary 
least squares (Isobe et al. 1990) used in previous studies 
(Schaefer 2007) is an alternative symmetrized approach 
(for other similar approaches see Kelly 2007). That ap- 
proach however ignores the measurement uncertainties 
during the fit which are taken into account in our sym- 
metrized maximum likelihood method. The advantage 
of expressing x^C'^i ^m) in terms of both the calibra- 
tion parameters a, b and the cosmological parameter 
Om is that it allows either fixing Qm and calibrating a, 
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b to be used for constraining other cosmological models 
(Schaefer 2007) or minimizing simultaneously with re- 
spect to all three parameters (Li et al. 2008b; Qi, Wang 
& Lu 2008). Even though the former approach may be 
useful as a consistency test it suffers from a circular- 
ity problem since it assumes a particular cosmological 
model to make the calibration and thus introduces a bias 
against alternative models. The later approach bypasses 
this circularity problem at the expense of increasing the 
uncertainties in the parameter determination since there 
are now three parameters to be fit instead of just two. 
Thus, in the next section we use the former approach 
(fixing Qn-i) to investigate the possible redshift depen- 
dence of the calibration parameters a, b but we use the 
later approach (variable ilm) when testing for the model 
dependence of the calibration. 



3 RESULTS OF CALIBRATION TESTS 

3.1 Redshift Dependence of Calibration 

In order to test if the correlation relations discussed in 
the previous section vary with redshift we separate the 
GRB samples in each case into four groups correspond- 
ing to the following redshift bins: z G [0,1], z £ [1,2], 
z G [2,3] and z G [3,7]. The number of GRBs corre- 
sponding to each correlation relation and each redshift 
bin is shown in Table 2. 

We now perform a maximum likelihood fit in each 
group after setting Sim = 0.26 obtained from the five 
years WMAP data (Komatsu et al. 2008), and deter- 
mine the best fit calibration parameters a, b with la 
errors and the quality of fit expressed through the min- 
imum (psr degree of freedom). We use a different 
value for the intrinsic dispersion aays for each type of 
correlation obtained by demanding a value of of order 
1. Thus we transfer the information about the quality 



Table 2. Number of GRBs in each redshift bin. 



Correl. Type 


z g [0, 1] 


2 e [1,2] 


z e [2, 3] 


2 e [3, 7] Total 


Epeak ~ Ej 


10 


8 


4 


5 27 


Epeak ^ L 


18 


15 


14 


17 64 


nag - L 


7 


13 


9 


9 38 


V-L 


14 


15 


9 


13 51 


TRT — L 


17 


15 


13 


17 62 



from to (Jsys (the smaller the required asya the better 
the fit). These results are shown in Table 3. We have ver- 
ified that the main features of our results do not depend 
on the choice of 0.-^ in the range Q.m £ [0.2, 0.3]. Despite 
of the fact that we use the maximum likelihood method 
instead of linear regression used in Schaefer (2007) our 
results for aj, bj in the full redshift range (last column 
of Table 3) are consistent at the la level with those of 
Schaefer (2007). 

The corresponding plots of the calibration param- 
eters aj, bj vs < z > (average redshift in each redshift 
bin) are shown in Fig. 1 for all five correlation relations. 
The slope of each aj{< z >), 6j(< z >) with la er- 
rors is also shown in Fig. 1. Notice that all slopes are 
consistent with zero at the 2a level. Thus, there is no 
statistically significant evidence for evolution with red- 
shift of the correlation relations. This result is to be con- 
trasted with the result of Li (2007a, 2007b) where the 
Amati et al. (2002) relation Epeak — E-y-iso was found 
to have statistically significant correlation with redshift 
for both parameters a and b. However the latter results 
are under dispute by the recent paper of Ghirlanda et 
al. (2008). Ifere we use an improved variant of this re- 
lation namely Epeak — E^ (Ghirlanda et al. 2004) and 
a somewhat different GRB dataset (Schaefer 2007) and 
we find no statistically significant evidence for correla- 
tion of a, b with redshift. It is also interesting to use 
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Table 3. Best fit parameters in each redshift bin 

Correlation Type 2 € [0, 1] 2 G [1, 2] 2 G [2, 3] z e [3, 7] Total 

Epeak-E^, 6 = 50.56 ±0.06 6 = 50.57 ±0.08 b = 50.63 ± 0.11 6 = 50.66 ±0.15 6 = 50.60 ± 0.04 

(o-sys = 0.15) a = 1.56 ±0.15 a =1.76 ±0.23 a = 0.84 ± 0.37 a = 1.45 ± 0.26 a = 1.56 ±0.11 



■^peak 



■ 0.30) 



b = 
a ■■ 



52.12 ±0.10 
= 1.55 ±0.14 
v2 = 1.11 



52.43 ±0.12 
= 1.99 ±0.33 
= 0.89 



b = 
a -- 



52.61 ± 0.27 
= 3.08 ± 1.19 
y2 = 0.40 



b = 
a - 



52.51 ±0.13 
= 1.14 ±0.41 
= 0.20 



b = 52.33 ± 0.06 
a = 1.75 ±0.11 
= 0.84 



Tlag 
{<Tsys - 



- L 
■■ 0.39) 



b : 

a - 



: 51.80 ±0.24 
-0.63 ± 0.31 
= 0.62 



52.32 ±0.14 
-0.94 ±0.23 
y2 = 0.70 



52.32 ± 0.18 
-0.93 ±0.25 
y2 = 0.63 



b -- 
a - 



: 52.52 ±0.19 
-0.31 ±0.26 
v2 = 0.13 



: 52.24 ± 0.09 
-0.87 ±0.12 
y2 = 0.87 



V- 



L 

■ 0.47) 



b = 
a - 



52.24 ±0.16 
= 0.96 ± 0.36 
X-2 = 1.22 



: 52.50 ±0.14 
= 2.34 ± 0.56 
y2 = 0.57 



6 = 
a - 



52.57 ±0.20 
= 1.50 ± 0.54 
y2 = 0.58 



6 = 
a - 



52.69 ±0.14 
= 0.03 ± 0.47 
v2 = 0.32 



: 52.45 ± 0.07 
= 1.30 ±0.21 
Y^ = 0.95 



trt 



-L 
■ 0.48) 



b-- 
a - 



■- 52.29 ±0.16 
-0.92 ±0.21 
v2 = 0.52 



b-- 
a - 



52.56 ±0.13 
-1.05 ±0.22 
y2 = 0.27 



: 52.43 ±0.13 
-0.81 ±0.26 
y2 = 0.56 



b- 
a - 



■ 52.59 ±0.11 
-0.43 ± 0.26 
Y^ = 0.19 



: 52.46 ± 0.06 
-0.92 ±0.11 
y2 = 0.87 



the results of Table 3 in order to compare the five cor- 
relation relations with respect to their quality of linear 
fit. Since the value of the intrinsic dispersion a^ys has 
been adjusted for each correlation so that — 1, the 
comparison can not be made by simply comparing the 
values of (per degree of freedom) in each case. In- 
stead we compare the required value of asys- Smaller 
asys corresponds to better quality of fit. According to 
this test, the correlation relations are ranked according 
to their quality as follows: 

(i) Epeak — E~/ {asys = 0.15), 

(ii) Epeak - L {asys = 0.30), 

(iii) Tiag - L {asys = 0.39), 

(iv) V-L {asys = 0.47), 



(v) TRT - L {c 



0.48). 



Clearly, the correlation relation Ep^ak—E^ provides sig- 
nificantly better fit compared to the other four relations. 



3.2 Model Dependence of Calibration 

In order to investigate the cosmological model depen- 
dence of the calibration parameters aj, bj we now al- 
low the cosmological parameter fim to vary in the 
X^(a, fe, Om) of eq. (5). In particular we minimize 
X^(a, 6, f^m) with respect to aj, bj for various nm € 
[0.2, 0.3]. We show Xmm(!^m) in Fig. 2 for all five corre- 
lation relations. As shown in Fig. 2 the best correlation 
relation {Ep^ak — E^) develops a minimum with respect 
to at fin, = O-SOtlfl (the la error is obtained by 
demanding < Xmm + 3.5S/N corresponding to max- 
imum likelihood minimization with three parameters). 
On the other hand, the for other four correlation 
relations are monotonic with respect to flm and seem 
to favor a flat matter dominated universe. We attribute 



this behavior to the large intrinsic scatter inherent in 
these correlation relations. The fact that these relations 
favor largo Sim may explain the fact that the best fit fim 
obtained by Schaefer (2007) using a combination of all 
five correlation was somewhat larger (fim ~ 0.39) than 
the corresponding values obtained using other more ro- 
bust cosmological tests. We thus argue that the mixing 
of the best quality correlation with the other four cor- 
relations plagued with large scatter may lead to mis- 
leading results and perhaps should be avoided. To this 
end, it is interesting to mention that systematic effects 
introduced by the flux limits of the GRB surveys (Lloyd 
et al. 2000) could potentially affect the above statisti- 
cal results. In order to check such a possibility, in the 
next section we perform a direct comparison between 
the GRBs with the Snia data. 



4 GRBS VERSUS SNIA 

An alternative model independent test of the GRBs as 
standard candles is their direct comparison with the 

best cosmological standard candles available namely 
SnIa. Thus we perform a cross-correlation analysis in 
the distance modulus space in the redshift range z € 
[0.17, 1.75]. The cross-correlation between two samples 
(in our case GRB and SnIa) is typically deflned by 
the following estimator (Peebles 1973; Efstathiou et al. 
1991): 

,Ngg 



$(m) = / 



Ngs 



1, 



(10) 



where Naa and Ngs is the number of GRB-GRB and 
GRB-SnIa pairs respectively in the interval [jj,, + Sji]. 
For a ^(/tt) = there is not any correlation between the 
two populations. Having known (from the x^— analysis) 
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Figure 3. Comparison of the distance modulus as a function 
of redshift, estimated by tiie GRB observations (solid points) 
and the a recent Snia dataset of Davis et al. (open points). 
Last Panel: The average cross-correlation function between 
the GRBs and SnIa as a function of the density parameter. 
The different points correspond to different luminosity indi- 
cators for the current GRB data and thus to different values 
of the fitted model parameters. Note that we don't plot all 
the error-bars in order to avoid confusion. 



the calibrations for every L or E~f, we can also derive 
[sec cqs (1) and (2)] the corresponding distance modu- 
lus jj,{z) = 51ogd_L(z) + 25. Note that Sfi is set at 0.75 
which corresponds to 10 number of bins. The robustness 
of our results was tested using different bins (spanning 
from 7 to 15) and we found very similar statistical re- 
sults. In the above relation / is the normalization fac- 
tor / = 2Ns/iNG - 1) where No and Ns are the total 
number of GRB and SnIa entries respe ctively. The un- 
certainty in ^{/j.) is estimated as = + /Ngs 
(Peebles 1973). To this end, we estimate the average 
cross-correlation function, which is given by: 



i=l 



(11) 



where Nb is the number of bins used {Nt — 10). 

In this paper we utilize the sample of 192 stan- 
dard candles (supernovae) of Davis ct al. (2007) in 
the 0.016 < z < 1.755. Due to the fact that the 
GRB observations extend beyond z > 0.17, we apply 
our statistical analysis to the following redshift interval 
0.170 < z < 1.755. In this case, the SnIa subsample 
contains 146 entries, Ep^ak — contains 17 entries, 
Epeak — L contains 32 entries, Tiag — L contains 19 en- 
tries, V — L contains 28 entries and trt — L contains 31 
entries. 

In Fig. 3 we compare the estimated GRB distance 
moduli (solid points) with those derived by the SnIa 
data (open points) as a function of redshift. Perform- 



ing, a standard Kolmogorov-Smirnov (KS) test to the 
distance moduli we find that the probability of consis- 
tency between GRBs and SnIa using the 

Epeak E^ 

relation is Vks ~ 0.4. Clearly, there is a strong indica- 
tion that the Epeak ~ E-, relation traces well the SnIa 
regime, a fact corroborated also by the cross-correlation 
test between GRBs and SnIa, which gives an average 
cross-correlation function of 0.38 ± 0.09 for fim = 0.30 
or a ~ 4(7 correlation signal. Notice that the cross- 
correlation analysis was performed for each value of Qm ■ 
Doing so, in the last panel of Fig. 3 we present, for each 
GRB subsample, the average cross-correlation function 
as a function of Om- It becomes clear that the best 
GBR tracer of the SnIa regime is the Epeak — E-, rela- 
tion. Interestingly, the cross correlation function peaks 
at f2m ~ 0.15 and flm — 0.30 respectively. However, 
due to small number statistics the Epeak — E^ test (17 
entries) provides much weaker constraints on fim than 
current SnIa data. On the other hand, the riag — L gives 
a relatively good correlation signal up to Om < 0.3 
although, the amplitude of the cross-correlation func- 
tion decreases rapidly as a function of Q.Ta than the 
Epeak — E^ case. Finally, the tbt — L, Epeak — L and 
V — L relations seem to fluctuate around zero (in agree- 
ment with the minimization results of Fig. 2). 



5 CONCLUSIONS 

We have investigated the robustness of the current GRB 
data calibrated as luminosity indicators with respect to 
two sources of biases: evolution of the calibration with 
redshift and dependence of the calibration on the as- 
sumed cosmological model. We have found no statisti- 
cally significant evidence for evolution of the calibration 
parameters a, b with redshift for any of the luminosity 
indicator correlations considered. However, our statisti- 
cal results do not exclude the possibility of correlation to 
be discovered in the future based on better GRB data. 
We have also found that the correlation Epeak — Ej has 
two important advantages over the other four correlar 
tions considered: 

• Its intrinsic scatter asys is less than half of the cor- 
responding scatter of the other four correlations 

• It can pick up the accelerating expansion of the uni- 
verse in a model independent way [it has a clear global 
minimum of x^{o-,b,ilin) at Om — 0.3] 

• It traces relatively well the SnIa regime. 

However, even the best GRB luminosity indicator cor- 
relation is currently not competitive with other cosmo- 
logical probes of the acceleration expansion since the 
cosmological parameter la errors (f2m = 0.30lJ 2") are 
more than an order of magnitude larger than the cor- 
responding errors obtained eg using SnIa standard can- 
dles and other geometrical probes (fim = 0.267lQ g^g). 
Therefore, even though the GRB data are currently not 
competitive with other cosmological probes of the accel- 
erating expansion of the universe this may well change 
in the future if the Epeak — E^ GRB dataset expands 
well beyond its current status consisting of only 27 dat- 
apoints. 
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